2210.09858v1l [cond-mat.stat-mech] 18 Oct 2022 


e 
e 


arXiv 


The vanishing of excess heat for 


nonequilibrium processes reaching zero ambient temperature 


Faezeh Khodabandehlou,! Christian Maes,! Irene Maes,” and Karel Netoény? 


‘Instituut voor Theoretische Fysica, KU Leuven, Belgium 
2Departement Wiskunde, KU Leuven, Belgium 


3 Institute of Physics, Czech Academy of Sciences, Prague, Czech Republic 


We present the mathematical ingredients for an extension of the Third Law of 
Thermodynamics (Nernst heat postulate) to nonequilibrium processes. The central 
quantity is the excess heat which measures the quasistatic addition to the steady 
dissipative power when a parameter in the dynamics is changed slowly. We prove 
for a class of driven Markov jump processes that it vanishes at zero environment 
temperature. As an application, the nonequilibrium heat capacity goes to zero with 
temperature as well. Main ingredients in the proof are the matrix-forest theorem 
for the relaxation behavior of the heat flux, and the matrix-tree theorem giving the 
low-temperature asymptotics of the stationary probability. From that arborification, 
the main new condition for the extended Third Law is dynamical and requires the 
absence of major (low-temperature induced) delays in the relaxation to the steady 


dissipative structure. 
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I. INTRODUCTION 


The Third Law of Thermodynamics was the culmination of a longer process of trying to 
understand thermal properties for macroscopic systems at low temperature. It is different 
in generality from the First and Second Laws, as it does not follow from the structure 
of equilibrium statistical mechanics. It is therefore also often called a Postulate, and one 
refers there to the version formulated by Walther Nernst around 1907, stating that the 


change of entropy in any isothermal process approaches zero as the temperature reaches 


absolute zero. Real progress in understanding came around the time of the first Solvay 
conference (1911), where the question of specific heats was intensively discussed and where 
it was suggested by Einstein that Planck’s theory of radiation should radically modify also 
the older kinetic-molecular theory”. In the end, and as the quantum revolution unfolded, 
measurements of heat capacity were finally understood as informing about energy spectra, 
fluctuations and occupation statistics, and the Third Law (in the later formulation of 
Planck) became a statement of the vanishing of the entropy of a perfect crystal of a pure 


substance. 


With the exception of some early thoughts by James Jeans (in his Dynamical Theory 
of Gases, 1904), nonequilibrium features modifying heat capacities have not been studied 
on a more fundamental level. In modern times, the nonequilibrium extension originates in 
constructions of steady state thermodynamics and in the notion of excess heat in particular, 
[214]. Nonequilibrium heat capacities were explicitly introduced and discussed in [5H8]. 
It was already observed there how their behavior as function of temperature may inform 
us about dynamical aspects as well. The main change from equilibrium is indeed that 
nonequilibrium heat capacities are able to pick up the dependence of the relaxation of 
heat and dissipated power on temperature and other parameters. The low-temperature 
behavior of heat capacities is therefore not only or no longer only informing us about static 
fluctuations but information about dynamical accessibility enters as well. That is seen 
most spectacularly in the behavior of heat capacities in the immediate neighborhood of 
absolute zero. As we will prove and illustrate in the present paper, we need a condition 
on the dynamics that ensures that there is little variation in the relaxation behavior of 
the dissipation. All states must remain “relatively well-connected” to the most probable 
zero-temperature state. In other words, not only is there a condition about absence 
of degeneracy of low-temperature stationary measures, but a dynamical accessibility or 


“no-delay” condition must be added. Or, even when there is a well-gapped unique invariant 


1 Its historical origin lies in the variational principle of Thomsen and Berthelot, which was an empirical 


precursor of the Gibbs variational principle, [I]. 
2 The problem of specific heats, the contradictions between the low-temperature empirical results and the 


equipartition theorem as e.g. applied in the law of Dulong-Petit, was indeed highlighted by people like 
Maxwell (already around 1860) and Kelvin, Jeans and Gibbs to name a few. 


measure at zero temperature, the heat capacity can still be different from zero (and even 
diverge). That happens at parameter values where the relaxation behavior of the heat gets 
pathological; we can speak here of a localization-phenomenon where there is a local delay 


in the relaxation to the stationary dissipation. 


Another difference with equilibrium, is that for nonequilibrium it remains relevant to 
consider small open systems. The temperature is not that of the system but of a thermal 
bath weakly coupled to the system. The system dissipates energy (heat) into that thermal 
bath, even when stationary. There is however an additional heat, an excess heat Q°° = Q, 
when considering the relaxation between two stationary nonequilibrium conditions. The 
main question of the present paper is to mathematically understand how the excess heat 
behaves when the temperature of the bath tends to zero. More precisely, we consider an 
averaged quantity of excess heat Q(T) along a process where parameters like temperature 
are changing, with I a curve in parameter space. For quasistatic changes in the parameters, 
we must prove that 

QT) = fea 
r 
where, depending on the variations, C, would correspond to a heat capacity or a latent 
heat. Important is that Q(T) is “geometric”, only to depend on the curve [ in parameter 
space and not on how fast temperature or other parameters change in time. Such things are 
known for currents in the quasistatic limit of infinitely slow processes; see [9]. Next comes 
the question how and why C) vanishes at absolute zero, to extend the Nernst postulate to 


nonequilibrium statistical mechanics. 


A physics presentation of the mathematical results that follow is found in [10]. In par- 
ticular, we motivate there our choice of excess heat and how our dynamics has quantum 
features. More illustrations and calculations of nonequilibrium heat capacities are found in 
[IHI3]. The present paper gives the rigorous version including all mathematical details. 

In Section [I] we start with the setup, definitions and a presentation of the main results. In 
particular, we specify the Markov jump processes with the physical interpretation of heat 
and dissipated power in a thermal bath, as follows from the condition of local detailed bal- 
ance, [I4]. From there follows the notion of quasipotential Vy (x) which is an expected excess 


heat when in state x as it depends on the parameters A. It appears in the quasistatic limit 


E | 0 of the expected excess heat Qe, 


QT) = lim Qe = fo (Va Vaa (1.1) 


where the (-), indicate the stationary measure at parameter values \ € R"*!. Our main 


result states that Q(A) vanishes as the temperature of the thermal bath tends to absolute 
Zero [Theorem [I.2]. Secondly, we consider the heat capacity C (8) = BDV à which is the 
variation of the excess heat with temperature, and we show that C(() tends to zero as well 
for inverse temperature 8 + oo [Theorem [1.3]. The conditions for both Theorems can be 
stated in graphical terms (such as from the matrix-forest theorem). Those conditions are 
illustrated in Section [I] with simple examples. 

The actual proofs of the main theorems start in Section [IV] It first concentrates on showing 
the geometric result from equation (I.1). The result is that the excess heat gets expressed 
as the path-integral (in parameter space) of V\dp\ where p$ is the stationary probability 
law. Section shows the boundedness of Vy in 8 f oo, and Section [VC] gives the 
low-temperature asymptotics of dp, —> 0. A summary of results and arguments is given in 


Section [V] The Appendix makes the explicit link with the matrix-forest theorem. 


II. SETUP AND MAIN RESULT 
A. Markov jump process 


Consider a simple, connected and finite graph G = (V, E) with vertex set V(G) and edge 
set E(G). Vertices are written as x, y,...; edges are denoted by e := {x,y} when unoriented 
and € := (x,y) is an oriented (or, directed) edge which starts in x and ends in y. 

A random walker with position X, € V at time t has transition rates k(x,y) > 0 for the 
jump x — y when {x,y} € E, and otherwise the rates are zero. By irreducibility, there 
is a unique stationary probability distribution p*(z) > 0,2 € V, solution of pL = 0, 
where L is the backward generator (|V| x |V|-matrix) having elements Ley = k(x, y) 
and Lge = = k(z,y). Stationary expectations with respect to p* are written as 
(PY = do. f(x) p*(x) for functions f(x), € V. We use the notation (f(X+)|Xo = x) for 


the conditional (process) expectation when starting in state x. 


We need to consider the dependence of the process on n+1 real parameters A = (3,a),a = 


(a;i = 1,...,n) with 8 > 0 interpreted as the inverse temperature of a thermal bath, and 


other parameters a € A for some open set A C R”. That leads to a notation with subscripts 


like in ky(x,y) and p$ that we use when relevant. We assume that all the transition rates 
ky(x,y) are smooth in À. Note here that we effectively deal with one single random walker 
on an arbitrary connected graph whose vertices represent possible many-body conditions. 
In that way, our work deals with general finite state-spaces with arbitrary interactions as 
long as they are smooth. 

For interpretation, we suppose that 


1 ky (a, y) 


B O ky (y, 2) = a(x, y) 


does not depend on 8 and q,(z,y) is obviously anti-symmetric. Following the physical 


condition of local detailed balance [14], qalx, y) can be interpreted as the heat to the thermal 


bath in the transition x — y. The expected instantaneous power gets then defined from 


Pro) = T ae y\ea(o59): (II.1) 


y 


By convexity (Py), > 0, and (P,), is called the stationary dissipated power. An important 


quantity will be the quasipotential Vy, a function on V defined as 
+00 
oe f TEA IEA E AN (11.2) 
0 


Given a smooth time-dependence X(t),0 < t <7, of the parameters, we call I its image 


as a curve in parameter space R* x A. For a quasistatic process, we write \° and consider 
a protocol where the system evolves under A*(t) := A(et),0 < t < e~'r on T. The e is the 
rate of change in the parameter protocol where À evolves. The limit where e | 0 will make 
the process to become quasistatic. For such a time-dependent process, at every moment t 
the distribution is pf, solving the time-dependent Master equation, 

Soi 


for the corresponding forward generator rit (transpose of Ly), L$ p5 = 0. 


= Len bi (11.3) 


For a given protocol A* at rate € > 0, the expected excess heat towards the thermal bath at 


inverse temperature ( is given by 


r/e 
Q: := f dt SO (Paen (2) pf la) — Pren(@) 4 e1)(2)). (11.4) 


The quasistatic asymptotics of (II.4) yields the “geometric” expression 


Proposition II.1. 
lim Qe = Q(T) = f dd (Vas), (11.5) 
in terms of the quasipotential (II.2). 


We call (II.5) geometric because the integral is over the curve T. It allows us, mathemat- 
ically and physically, to define the heat capacity 


C(B) = 6 (Ns), (II.6) 


which fixes all parameters œ in the variation. 


B. Graph elements 


We recall some standard notions from graph theory. Remember that we have a connected 
and finite graph G = (V,€) with vertices x, y,... and edges denoted by e := {x,y} when 
unoriented and € := (x,y) if oriented (or, directed) from x to y. 

A subgraph H of G is a graph for which V(H) C V(G) and €(H) C €(G). We call it a 
spanning subgraph if V(H) = V(G). An oriented subgraph, meaning a subgraph where 
every edge has a direction, is shown by H. 

A path on the graph G is an alternating sequence of vertices and directed edges where all 
visited vertices are different. For example, (x, (x,y), y, (y,z),z) is a path from x to z. An 
oriented loop denoted by f is an alternating sequence of vertices and directed edges where the 
initial and the final vertex are the same and the other vertices are all different. For example, 
(x, (x,y), y, (y, z), 2, (2,2), x) is an oriented loop. A loop denoted by £ can be obtained from 
an oriented loop by removing orientations of edges. 

A rooted tree is a connected oriented subgraph of G without any loop. In a rooted tree there 
is always a special vertex x such that for every vertex y, there is a path from y to x and 
there is no edge going out from x. We call this special vertex x the root of the tree. We use 
the notation 7, to denote a spanning tree rooted in x. A tree in a graph G is a connected 
subgraph of G which can be obtained from a rooted tree in G by removing all directions of 
the edges and only have unoriented edges. We use the notation 7 to denote a non-oriented 


spanning tree. 


Since tree-loops and tree-loop-trees are especially important in what follows, we introduce 


them in the following definitions. 


Definition II.1. (a) An oriented tree-loop is a connected oriented graph consisting of an 


oriented loop and rooted trees such that the root of every tree is located on the loop; 


see Fig. 


(a) counter clockwise (b) clockwise 


FIG. 1: Two oriented tree-loops which only differ in the orientation on the loop. Remark 


that in every tree-loop the roots of the trees are located on the loop. 


(b) A tree-loop is a connected (non-oriented) graph that can be obtained from an oriented 
tree-loop by removing the orientations of the edges; see Fig.|2} We denote the set of all 
unoriented spanning tree-loops in the graph G by H, the set of all spanning tree-loops 
such that the edge e is part of a tree by He and O(He) denotes the set of all oriented 


tree-loops including the edge e located on a tree. 


FIG. 2: A tree-loop. Any number of trees may be connected to the loop. 


Definition II.2. (a) A tree-loop-tree is a graph consisting of two disconnected parts, one 


of which is a tree and one of which is a tree-loop. We can make it by removing an edge 


from the tree-parts in a tree-loop; see Fig. The set of all spanning tree-loop-trees 


made by removing the edge e from the tree-loops in the set He is denoted by H™. 


FIG. 3: A tree-loop-tree made by removing the edge e from a tree-loop. 


(b) An oriented tree-loop-tree is a tree-loop-tree where every edge has a direction. The 
directions of the edges are such that the tree-loop part is an oriented tree-loop and 
the tree part is a rooted tree; see Fig. If H € He is a spanning tree-loop, then 
O(H)) denotes the set of all oriented spanning tree-loop-trees which are made by 
removing the edge e of H and giving all possible orientations (there are two different 
orientations possible in the loop and different orientations possible in the separate tree 
depending on the location of the root in that tree) to the resulting tree-loop-tree. We 
define o4(€) = —1 if € is oriented towards the loop of H and og(€) = +1 otherwise. 


In this work an (oriented) tree-loop-tree means a spanning (oriented) tree-loop-tree. 


FIG. 4: An oriented tree-loop-tree made by removing the edge e from a tree-loop and adding 


orientations. The orientation in the loop is counter-clockwise and the tree part is rooted in r. 


We define the weight w(€) for an oriented edge as w(€) = k(€) . The weight of an oriented 
subgraph H with €(H) # 0 is 
w(H)= || wv. (II.7) 


If there is no edge in H but it still has vertices, then we define w(H) = 1. The weight of 
the empty set is w(0) = 0. We define W = 90.4, 357 w(Tz) as the sum over the weights of 


all rooted spanning trees in the graph. 


C. Main results 


To state the main results, we introduce 


p(z, y) = Jim, 5log k(x, y) 
and for an oriented subgraph G of G, write (9) := 0. neeg) 92:2’). Define 
(2) = maxo(T), 6" = max max d(T.) = 4(0"). (11.8) 
where we call x* a dominant state. Here is the first condition. It consists of two subcondi- 
tions: 


Condition 1a: For every edge e in the graph G and for every element H of the set O(H) 
(see Definition. 
oH) < ¢*. (T1.9) 

That condition will be interpreted as a ‘no-delay’ condition; see Example 
Condition 1b: There is a unique z* € V independent of a € A so that with \ = (8, a), as 
P To, 

AlL) = za (1 + 0(1)) (11.10) 
where the correction o(1) is going to zero uniformly in a € A. 


The main result is a Nernst heat theorem for the given context. Recall Proposition (IT. 7] 


Theorem II.2. Recall equation (11.4). Under Condition 1, along any curve T in A (at fixed 


B), 


lim lim Q; = lim Q(T) = 0. Thi 
Ae eh Oe FO) eon 


In the next Section, it is shown how Condition 1 is natural, and physically meaningful. 


We next consider the heat capacity (II.6) in the limit 6 + oo. We fix the other parameters 
a € A. We need the second condition: 


Condition 2: limg;. p*(x) = p$ (x) such that 


p(z) = p(x) (1 + O(8*)) (11.12) 
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with O(8~*) for 8 f oo uniformly in a. 


Theorem II.3. Recall equation (IL6). Under Conditions 1a and 2, 


lim C(B) = 0. (11.13) 


For Condition 1b and Condition 2, it is useful to remember the Kirchhoff formula 
for the stationary distribution p, 


w(2) 


p(t) = -y w(x) = > wT) (11.14) 


where W is the normalization and the sum in the weights w(x) is over all spanning trees T 
in the graph. The weight w(7,) is the product of weights (II.7) over the oriented edges in 
the spanning tree 7 when orienting the tree towards root x. As shown in [I6], this formula 


leads to a low-temperature asymptotics, [17] [8] 
p3(x) = Pa + Ole) (11.15) 


for a ô > 0 uniformly in a € A. Therefore, Condition 1b is implied when the maximizer 
of ¢* in equation (II.8) is unique but for Condition 1b the correction does not need to be 
exponentially small in large 8. For Condition 2, we do not even need nondegeneracy of the 


limiting stationary measure. 


We end the section by mentioning the key to the main results which is a graphical 


representation of the difference-quasipotential over any edge € = (x, y). 


Proposition II.4. The difference-quasipotential V (€) = V(x) — Va (y) (from (II.2})) can be 
decomposed as V (€) = Viree(€) + Vioopl€) with 
1 
Viree(®) = a= 2o arhe >y) D WT) (11.16) 
T ueV 
where the first sum is over all spanning trees, q(x — y) is the sum of q(z, 2’) over the path 
from x to y located on the tree T, and 
1 — — 
Vioow(®) = yr 2 aO J uE) (11.17) 
HEHe HEO(H(®)) 
in which the oriented loop £ is the one in H and oy(€) = +1 depending on whether € looks 


away or towards the loop (see Definition b)). 
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Remark: There are other representations of Viree(€) such as from observing that 


Viree(€) = g(€) + 5 > (ar(a > y) + a(y, 2) NWT) (11.18) 
T 


ueV 
where the second term makes a loop. That is useful for interpretation as it reduces to 
V(€) = Vrele) = g(€) in the case of an equilibrium dynamics. More generally, only 


heat-carrying loops contribute. 


Proofs are collected in Section We first give some physics explanations and illustra- 


tions in the next section. 


II. EXAMPLES AND ILLUSTRATIONS 


This section is meant to clarify the graphical conditions of the main Theorems |IT.2HI1.3 
It is not essential for the proofs and it can be skipped at first reading. We first illustrate 


the notation with some simple examples. 


Example I1.1. Consider the graph G in Fig. 


x y 


u Z 


FIG. 5: Graph G made by three loops. 


The transition rates are 
k(y, x) = k(y,z) = k(x,u) = k(u, z) = e’ 
and the rates over all other edges are equal to one. Definition [I.8] gives 


oy, x) = oly, z) T olz, u) = olu, z) =-], 
olz, y) = b(z,y) a o(u, x) = b(z, u) = (z, x) B olz, z) = 
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To find ¢@* we need to look at all rooted spanning trees in the graph. The spanning trees 
x Y x y 1] a 
u Z uU Z u Z u Z 
x y X y NS : | 
u Z u Z u Z u Z 


FIG. 6: Spanning trees of the graph in Fig. b See for more details. 


are in Fig. [6] 


We find 


and thus ¢* = 0. Denote 


e1 := {x,y}, CS {y,z}, e3:= {z,u}, e4:= {x,u}, e5:= {x,z}. 


To check Condition la we also need to find all spanning tree-loops in the graph which are 


shown in Fig. 


u Z u Z u Z u Z u Z 
Hı H3 H3 H4 Hs 


FIG. 7: All tree-loops in the graph of Fig. b 


For every edge, we need to look at the set of all spanning tree-loops that include this 


edge on a tree, 
He = {H}, Hez = {H2}, Hez = {H3}, He = {H4}, Hes = 0. 


The next step is to construct for every possible edge the set of all oriented tree-loop-trees 


O(H) (but that is not possible for the edge e;). For example, if we look at the edge e1, 
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we need the set O(H!*), which consists of two elements depending on the two possible 


orientations in the loop of H; see Fig. 


x ev x av 


FIG. 8: The possible oriented tree-loop-trees created by removing the edge e, from H4. 


Therefore, O(H®®) = {Hi1,H12}, where (Hi) = 0 and ¢(H12) = —2. We now do 
the same for the other edges and get, 


O( HS) = {H21, H22}, ¢(H21) =0, d(Ho2) = -2 
Oa) = {H31,H32}, ¢(H31) =—1, ¢(H32) = —1 
O( Hi) = {H41,H42}, (Haa) =—1, 6(A42) = —1. 


We see that Condition la is satisfied. 
To check Condition 1b, we can use the Kirchhoff formula which expresses the stationary 
distribution in terms of weights of rooted spanning trees; see e.g. [18]. Here, the state y is 


the unique dominant state: x* = y, unique maximizer of (z) in (II.8). 


Next we give a counterexample to Condition la, which can be interpreted as a “no-delay 
condition”; for details of such a physical interpretation see the earlier discussion in [$| and 


the recent [10]. 


Example 11.2. Consider the graph in Fig. [9|made by a centered triangle such that each state 


is symmetrically connected to the state in the center. 


= y 


FIG. 9: The graph has four loops. 


14 
A random walker is moving on the graph with transition rates 
k(x,y) = k(y, z) =k(z,2) =v, k(y,x) = k(a,z) = k(z,y) = ve * 
k(x, u) = ku) =k(z,u) = pe, k(u,2) = k(u,y) = k(u,z) = we OH? 


where £, a, A > 0. The spanning trees with the largest weight are rooted in u; they start 
from a state on the triangle and visit the next state on the outer triangle in clockwise 
direction before going to the center. Therefore, @* = —a. 


There are 12 spanning tree-loops, shown in Fig. 


x x 


x x 
VA y Z y Z Y Z 2 y 
Hı Hp H3 Ay 
x x x x 
u u 
á y 2 yY z t% z "~y 
He Hy Hg H; 
x x x x 
u PAN /N 
Z y z y Z yY 2 y 
Ho Hio Hy, Hy2 


FIG. 10: All spanning tree-loops in the graph of Fig. b 


Putting 


e1 := {x,y}, eH {y,z}, e3:={z,£}, e4:= {x,u}, es:= {y,u}, ee {2,u} 


we have 


He = {H3, He}, Hez = {H4, H7}, Hes = {H,, Ho}, 
He, = {H2, Hi2}, Hee = {H;, Hio}, Hes = {Hs, Hi}. 
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The @’s for all small loops in the clockwise direction are equal to —(2 a + A) and in counter- 
clockwise direction equal to —(2a + A +€). The set 1‘) contains two tree-loop-trees made 
by removing the edge {x,y} from Hz and Hg. If we give orientations to both, we get four 
oriented tree-loop-trees for which the ¢ value is either equal to —(2a+A) or to —(2a+A+e). 
That scenario repeats itself for H) and H). Next, we take an edge connecting the outer 
triangle to the center and look at the sets H'4) = H) = Hle) all containing two elements. 
Giving orientations to these two tree-loop-trees we get four oriented tree-loop-trees. The @ 
of the large loop in clockwise direction is zero and in counter-clockwise direction —3¢e. We 
see that Condition la is satisfied for the edges e1, €2,e3 but not for the edges e4,e5, eg. In 
fact, it is easy to convince oneself that V(x), V (y) and V(z) are diverging for a > 0 while 
V(u) is uniformly bounded only when A > a. That is interesting because it shows that the 
conditions of our Theorems are not necessary: for a < A the heat capacity and the excess 
heat do go to zero at absolute zero. The reason is that the divergence of the quasipotential 
V(x) = V(y) = V(z) in £ fT œ is slower than how their stationary probabilities go to zero. 
For a > A also V(u) is diverging and the heat capacity as well. That again shows the 
physical content of the condition la: if a > A, there is a high barrier between the states 
x,y,z on the one hand and z on the other hand. Starting say from state x shows much delay 


in reaching the dominant state u. 


IV. PROOFS 
A. Quasistatic analysis 


Note that always >>, L'p(x) = X, L1(x) p(x) = 0 because L1(xz) = 0. Hence, when 
Lip = v, then >>,,v(x) = 0 must be true necessarily. It means that the (nonexisting) 


inverse 1/L' can only be defined, if at all, on those v with }>,. v(x) = 0. It implies 


D dt e™[f(x) — (f)*] v(x) = D dt e"! f(x) v(x). 


Therefore, we can actually define (1/L') v via the requirement that for all functions f, 


LIOGO- h ae - Wo) 


x 
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z A dt [ProbLX; = 2|Xp = y] — pE] v(y) 


which is indeed well-defined (and we do not bother introducing new notation for that 


restricted inverse). 


We want to solve the time-dependent master equation perturbatively, 
Pi = Ple HE p +0(e°) 


as everything is smooth around ¢ = 0. Since time gets rescaled by € we know automatically 
that 2 F® = O(e). On the other hand, 


3 pe) 2 


a BE lo; — Pye)| =L 


ea =E À(et) t VAPA) 


Hence, 


. 1 
FO = Net)» Va. 
L 
AÀ 
As a consequence, 


Tle 1 
f (Pi — P\et) dt = [aX Sve +0 (IV.1) 
0 T Ly 
where the integral on the right-hand side is purely geometrical, invariant under a 


reparametrization of time t > 7(t) for a smooth function y with y(t) > 0. 


Suppose we have a function f, on V. We can always write fa = —LaVa + (fy)}. Define 


the excess f$*° := fy — (fy)5, with time-integral 


em fo a (x 


In the case of f, = Py, from (I.4) we have Qe = (ey. 


Proof of Proposition |[IT.1} From (IV.1), 
T/E 


lim (de) = Tim i dt N (fren (x)p ile) — faen (2) 0 (et) (2)) 


fo TA anne 
= - far. Late z) Vap (z (IV.2) 


= | dà- vnose -2 P (£)V Val (IV.3) 
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and 


lim ( EA Tt jat) = (Vaa = (Vido + fe i (Va Vaa (IV.4) 


when a, b are eee the initial and final points in the protocol. But, (Vs), = (Va)a = 0 


so that we end up with 


tim ( f ODE IEAA (IV.5) 


eL0 
which indeed means that the left-hand side is geometrical and in the case fy = Py, that 


identifies the thermal-response coefficient 


Cy = (VaVa), (IV.6) 


as wished. (Some more details on the above are given in [8]). 


The equality (IV.2) shows that the excess heat under a quasistatic protocol equals 


Tle ë 
im( | ety =- fa: LORO 
eļ0 0 
We have thus reached 
lim Qe = i a: DOV) Voile (IV.7) 


As a consequence, to prove Theorem [IT.2| 2| the anes boundedness of the quasipotential V 
combined with the vanishing of V)p}\(x) as 6 + œo suffice. We start with the boundedness 


of V in the next section. 


B. Boundedness of the quasipotential 


Proposition gives a graphical representation of the difference-quasipotential V (e). 
The main inspiration and tool is the matrix-forest theorem. The proof is in Appendix [A] 
We start however with a direct proof that the difference-quasipotential is indeed derivable 


from a potential. 
Lemma IV.1. For all loops l in the graph G, and all oriented loops £, Y` z-z V (E) = 0. 


Proof. Let @ be the oriented loop ((21, £2), (£2, £3), . -. , (En-1, Zn), (£n, 21)). We start by 
showing that J scz Viree(€) = 0. Put &n41 = 21. Then 


Viree(€ = => Ss" w( w(Tu 103 aT (Xi = zu) . 


Ect T uev i=1 
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Since for every fixed J and every three vertices x, y, z, qr(£ > y)+qr(y > z2) = qr(z > z), 
the sum > 77) @7 (a; > Ti+1) = 0 and thus } zez Virec(€) = 0. 

Consider a graph G which has only one loop ¢. The edges located on the loop in the 
underlining graph are not participating in some tree of tree-loops of the graph G. So then 
He = Ý for all € € £ and then ))¢7 Vicop(@) = 0 (see (II-17)). 

Let us assume that the graph has more than one loop. In that case, there are edges on 
a loop for which there is a chance to participate on the tree part of a spanning tree-loop. 
We fix an oriented loop @ and sum over its edges, notice that there is more than one path 
between every two states on the loop /. The oriented loop @ has at least three edges. We 
call two of them € = (z,x),@ = (x,y). For the other loop there is a spanning tree-loop 
He, % e2 such that e, is located on a tree and g; is leaving the loop so that om., = +1. 
Consider a tree-loop made by removing the edge e; from He, and adding the edge e2. The 
new spanning tree-loop is denoted by He, % e;. When € in He, is oriented towards the 
loop, then oy(@,) = —1, while H = H2), That scenario repeats for every edge in £: we 


claim that for every H‘ with the sign o(@) there is a same HC” with the sign —oz(e’) 


such that the tree-loop-trees are the same and therefore X cT Vioop(e) = 0. 
Next, we give a graphical representation of (II.1) in the following 


Lemma IV.2. The average of P(x) in is 
1 a 
(P) = w X wa (IV.8) 


HEO(H) 


where O(H) is the set of all oriented tree-loops in graph G. The loop £ is same as the loop 


in H. 


Proof. See Appendix [B] 


The main step for proving Proposition is the next proposition. 


Proposition IV.3. For every x € V, 
XO ke) (-V(@) +4@) = (P). (IV.9) 
é=(z,-) 


Proof. We refer to the decomposition (II.16)) and (II.17) of V(é). We start by calculating 
dex(x,) K(@)Viree(€). If the edge e = {x, z} is located on a spanning tree T, then gr(x > 


e 
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z) = q(z,z), otherwise gr(x —> z) = (qr(x —> z) + q(z,2)) + a(x, z), where the first two 
terms give q over an oriented loop made by the path from x to z on T together with the 


edge (z,x). It follows that 


HWM E HOLY w(t) a0) 
é=(z,-) s 
+y : (> SY w(T, ) (qr(x = z) + q(z,x£)). 


e=(z,-) YEV e€T 


In the first line we can use that Dorey Xor w(T,) = W and the first sum is thus equal to 
Že-(,) K(@)a(@). Now look at the second line. Take € = (x,z) and Jy. We look at two 
cases: y= x andy # x. If y = x, then by adding the edge € to 7, an oriented spanning 
tree-loop is made (but q over the loop is in the opposite direction of the oriented spanning 
tree-loop). Every oriented spanning tree-loop with x on the loop can be made in this way. 

If y # x and the edge (x, z) is not in 7,, there exists another edge such (x, 2’) in 7y. In the 
case that y is on the path from z to z in the tree, by adding the edge € = (x, z) to Ty, again 
a loop is made. If (x, z’) is on this loop, then there exists another spanning tree J’ made by 


removing the edge {x, z’} from T and adding the edge {x, z} such that 
k(x, z)w(Ty) (@r(@ > z) + 9(z,2)) + k(x, z JwT) (ar: (£ > 2’) + g(2’,2)) = 0 


where k(x, z)w(Ty) = k(x, 2’)w(T,). We thus have 


do FCVireol@) = X kOe) — = = > 2i 


@=(c,-) @=(c,-) W aa HeO(H 
1 
-W DA k(é A Su (qr(z > x) + q(x, z)). (IV.10) 
e=(x,) P e€T 


Now we look at the loop term Vioop(€). Let € = (#,z). Take an arbitrary H € He. If 
x is not on the loop part and not the root of the tree part in H, then x has at least two 


neighbours. Let e} = {x, 21} and eg = {z, z2}, then 


k(t.) XO on(@)q(Qw(H) +k(a,22) X, on(@)q(2)w(H) = 0 


HEO(H1)) HeEO(H2)) 
x is no root x is no root 


because oy has different signs. For more than two neighbours on the tree the result is the 


same. In case that x is the root of the separate tree of some H € O(H), then by adding 
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the edge Z to H, an oriented spanning tree-loop is made. We thus have 
1 
LOM) =p D 


HeHe HEO(H 
he HEO(H) 


oy DOD DL on) a w(F). (IV.11) 
E- 


HEHe HeEO(H)) 
xel 


In the second line, k(@)w(H) (see Fig. |11a) is the same as k(e’)w(7,) (see Fig.|11b), where 


€ = (a, z') with z’ such that e’ is on the oriented loop part of H and r is the root of the 
tree part of H. The second line of (IV.11) thus corresponds with the second line of (IV.10) 
(see Fig. and as a result we get 


` k(é) V(e) = DH k(é Vioop(€ e) + Viree(€ )) 
é=(z,-) 
= Y k@)q(@) - D9 dw. 
e=(z,-) H HeO(#H) 


The proof is finished by using Lemma |IV.2] 


mee 


FIG. 11: The second line of (IV.11) corresponds with the second line of (IV.10). 


(b) CUT 


Proposition IV.4. Under Condition 1a, the quasipotential V (€) is bounded for every edge 


é, uniformly for B T oo. 
Proof. Look at the tree term (I1.16) in Lemma (II.4). The terms qr(x — y) are bounded 


and W > w(T,) for all y. Therefore, Viree(€) is uniformly bounded in temperature. 


Consider next the loop term (I.17) in Lemma (I.4). There, H € O(H) is an oriented 


tree-loop-tree created by removing the edge e from the spanning tree-loop H and then giving 
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direction to the edges. In the low-temperature asymptotics, the weights are 


w(H) ~ ef} 


and 


W x~ ` ePP(Ty) w eho. 
Ty 


On the other hand, q(£) is bounded. Hence, 


is bounded if (H) < ¢*. 


Now we can finish the proof of the uniform boundedness of the quasipotential. 


Proposition IV.5. Suppose Condition 1a. The quasipotential Vy = V of (II.2) is uniformly 
bounded in p + co. 


Proof. If there is a unique state x* as in Condition 1b which gets all probability and for 
which still V(x*) is diverging, then there must be an oriented edge € where V (e) is diverging. 
The reason is as follows. 
Consider the case where V(az*) is diverging at zero temperature. We know that for all 
temperatures 

0 = (V) =V(a") p°(2") + Do[V(2") — D> vE). (IV.12) 

zp e:a*—z 

Now consider the limit 6 ¢ co and assume that all V(é) remain bounded. We also have 
p*(x*) + 1, p*(z) — 0. That is not compatible with and hence if V (x*) is diverging, 
then V(e) must be diverging as well for some e. That is however impossible under Condi- 


tion la. So then, the boundedness of the quasipotential over all oriented edges gives the 


boundedness of the quasipotential in every state. 


To prove Theorem [IT.2] it remains to show that Vp }(x) is vanishing as 6 f oo. We do 


this in the next section. 
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C. Low-temperature stationary measure 


As for equilibrium, the behavior of the static fluctuations also matter greatly at low 
temperature. Since we are dealing with finite systems (even small, as is relevant for station- 
ary nonequilibria), and no estimates on the behavior of the excess heat as function of the 
number of particles are attempted, we only need to worry here about the low-temperature 
asymptotics of the stationary probability. We refer to for a systematic study of 
the low-temperature structure of stationary measures for Markov jump processes; see also 


T. 


Recall that the parameters a are unrelated to temperature; see Section 
Proposition IV.6. Assume Condition 1b in (II.10). Then, limg Vap\(x) = 0 for alla € V. 


Proof. We have assumed smoothness in a@ for every 6 and the number of vertices is finite 
with stationary distribution p$ for an irreducible Markov process. Since x* in (II.10) is 
independent of a € A, and the correction o(1) is going to zero uniformly in a as well, 


we can exchange the limits “limg Va = Vq lima’ 


p Î o. 


and have the required Vap$(x) > 0 as 


Together with the uniform boundedness of the quasipotental V as shown in Proposition 


TV.5| that concludes the proof of Theorem [IT.2] 
Proposition IV.7. Assuming Condition 2, it follows that limg P SpE) = 0 forallx E€ V. 
Proof. We must understand 
halz) = tim fe (x) LEV 
a p 6 d8 PA ’ . 
Condition 2 assumes that there is a function calx), x € V, so that for all sufficiently large £, 


ex(2) — ca(x)| < O(B™) 


so that limg p\ = Ca, and dividing by 1/8 still goes to zero, so to have h(x) = 0. 


From the two previous sections we are done with the proof of Theorem as well. 
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V. SUMMARY AND CONCLUSION 


We have considered nonequilibrium Markov jump processes on a finite connected graph 
under the interpretation of local detailed balance. The main object is the excess heat Qe 
which is computed in the quasistatic limit £ | 0. It gives rise to an expression in terms of 


the quasipotential V} defined in (H.2): 
lim Qe = fe NOVa) Vapi (z). 


where pġ is the unique stationary distribution. 

Our main result is the formulation of sufficient conditions for that quasistatic excess heat 
to vanish as the temperature goes to zero. First, the low-temperature asymptotics of the 
stationary distribution ensures that V\p$(x) — 0. Secondly, we need the boundedness of 
the quasipotential V), uniformly in 6 f oo. It is in fact easier to look at the differences of 
V, over an edge, but the main ingredient is the matrix-forest theorem, which we further 
discuss in the Appendix [A] It gives a graphical representation of the quasipotential, and it 
allows the evaluation of the zero-temperature limit. In particular, we need to focus there on 
a decomposition of V, and the contribution of loops. The result can be inferred by explicit 
inspection, where the condition arises more specifically from the analysis of tree-loop and 
tree-loop-tree contributions at low temperatures. For further physics interpretation we refer 
to [10]. 


As a conclusion, a nonequilibrium version of the Nernst postulate is obtained. 


Acknowledgment: KN thanks A. Lazarescu and W. O’Kelly de Galway for previous in- 
spiring discussions on the subject. The work was concluded while authors FK and IM visited 


KN at the Institute of physics in Prague. They are grateful for the hospitality there. 
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Appendix A: Matrix-forest theorem 


As it may seem less clear how (IJ.16) and (1.17) arise, we give here their origin from the 


matrix-forest theorem. 


Consider a continuous-time irreducible Markov process on a finite state space V charac- 


terized by transition rates k(x,y) > 0 for x,y E€ V. Let L be the backward generator, 
Ly =k(z,y), cAyand z,y €V 


Consider N 
Vale) = | KIX = 2) — (oat (A.2) 
where we subtract the asymptotic stationary value 
(glz) | eo = 2) “> (9). 


Putting f(x) = g(x) — (g), 


Ve (ae) = Vela) = f * l(a) dt (A.3) 
so then 
, —]1 
Va) = lim oyn) 


where — is the resolvent-inverse of the backward generator L (see[13}). 

For our purposes, the set V is the vertex set of a connected graph G, and the transitions 
happen over its edges. We use to obtain a graphical representation of V+; put 
V := Vy. A spanning forest is a collection of trees that forms a spanning subgraph. Define 
the set Fp” to be the set of all spanning forests in G with k edges having the properties: 
every tree in the forest is a rooted tree, y is the root of one of the trees and x and y are in 
the same tree (so there is a path from x to y). Define Fp as the union of sets F°” in graph 
G. 

Proposition A.1. 
w (Fay) 


V(x) = ` uoa 
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Proof. From [20], for all a > 0, 


(>) Se 


where w(F, “”) is the weight of the set Fg” 


fo) = 3. I] #2’. 


Fer, Y (z,z')EF 


Thus, 


Deo o a A ) 
= lim 
2 a—oo kw(Fr) 
E C EEEN OEN og OF) 
2 dia yo ` a kw (Fy, + 2 im a atw( Fp) ` 


Next, use that (f) = 0 and hence, `, p(y) f(y) = 0. For the stationary distribution p we 
use the Kirchhoff formula (11.14). Therefore, >>, w(F,\’) f(y) = 0, and the second term in 


n 


the last line above is equal to zero. To continue the calculation, 


De Tan (Fr) Fy) 
lim 
fo Da a ow w( Fp) 


Define F*~Y := F?TY as the set of all spanning forests consisting of two trees. Remark 
that F,,_, is the set of all rooted spanning trees and W is the sum over the weights of all 


rooted spanning trees. 


Corollary A.2. A solution V with (V) =0 of LV = —f with (f) =0 is given by 


Ve) = zp ele) Fl). (Ad) 
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1. Quasipotential for a specific source 


Here we consider the quasipotential of (A.4) for a specific source. We focus on the case 


f(y) = Py) — (P) in (AA), where Ply) = E, kly, x)a(y, x) and q(2,y) = —q(y, x). We 
define 


doy WFP) Ply) ae) 
Pry. SY (Pes y 
V” (x): wW , Vasa (P) wW ; 
So then, the quasipotential can be written as 
V(x) = V? (x) — VP (x) (A.5) 


Proposition |A.6| shows that the quasipotential in (A.5) can be decomposed into two terms, 
one related to spanning trees only and the other containing loops. 


FY denotes a forest in the set of F°7¥. Write ky. := k(y, x) and q(y, £) := qyx. 


Lemma A.3. For each x on the graph G, 


Samu, SWF) e (A.6) 
y zy PREF TTY 
(zy) FY 
Proof. The product w(F™) ky, is the weight when adding (y,z) to the forest F®™. Let us 
consider forests FY and F'** which have different directions for the edge {z, y}: (z,y) is in 


the forest F” € F?” and the edge (y, z) is in the forest F'** € F”?7 . Adding the edge 
(y, z) to the forest FY is the same as adding the edge (z, y) to the forest F7; see Fig. 


(a) FY U (y, 2) (b) F**U (z,y) 


FIG. 12: Adding the edge (y, z) to F** and adding the edge (y, z) to the forest F”. 


and 


w(F™) k(y, z)a(y, z) + w(F*) k(z, y)alz, y) = 0. 
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So then 


> uF ryP iy) = > wl) SO ene 
= > y w( F”) kyz qyz. 


zy F?YEF®™Y 
(z,y)ERCY 


Lemma A.4. By adding the edge (y,z) E€ G to the forest F” € F*Y where (z,y) ¢ F”, 


the new graph is either a rooted spanning tree or an oriented tree-loop-tree. 


Proof. Consider the forest F”. It has two trees: one is rooted in y and the other tree is 
rooted in some vertex r. Vertices x and y are located in a same tree and the edge (z, y) is 
not in the forest. If the vertex z is on the same tree with y, adding the edge (y, z) creates 


an oriented tree-loop-tree, Fig.|13(a), and if z is in another tree then a rooted spanning tree 


is created, Fig. [13{b). 


(a) F™ U (y, 2) (b) F= U (y, z) 


FIG. 13: The edge (y, z) is added to the forest F”. 


We need extra notation for the next Lemma. Put K, for the set of all tree-loop-trees such 
that x is located on the tree-loop part. If x € Kz, then O(x) is the set of oriented tree-loop- 
trees for all possible directions in x and = denotes an element from the set O(K). £ is the 
set of all non-oriented loops in the graph G and Ky, denotes the set of all tree-loop-trees 


including the non-oriented loop £ and zx in the tree-loop part. 


Lemma A.5. VP (zx) splits into two parts; one where we only sum over spanning trees and 


one where loops are present: 


V? (x) = Viree(®) + Vicop(2): (A.7) 
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Here, the first part is equal to 


j.i 


Viree(®) = 37 X S w(%) are >y) 
yEN T 
and 
qr(@ > y) := q(u, u’) 
(uyu!)e (ay) 


where (x G y) is the path from x to y in the spanning tree T. The second part is equal to 


eE E E NA 


LEL KEKi x REO(K) 
where € is same as the loop in R and q(€) = ued alu, u’). 


Proof. Consider the first case in Lemma [A4] there are different possible tree-loop-trees for 
every non-oriented loop in the graph G. For every edge (u, u’) located on the loop, there 
are two paths between u and u’. Removing the edge (u, u’) from the tree-loop-tree gives a 
forest where u and wu’ are in the same connected component. There is a forest F” such that 


uw € F™ and (u, u’) ¢ F™, and 


(u,u’ Ee 
Let us go back to Lemma |[A.4] For a rooted spanning tree 7, in graph G, there is a unique 


path (x = y) from z to y on this tree. Then, 


w(Ty)ar(@ > y) = w( Fe) Kav tuw (A.8) 
(uu!) Eley) 
where F% € F**". We have used the fact that by omitting an edge {u,u’} located on the 
path (a a y) in tree J,, a forest including two trees, T, (toward u) and T}, is created. 
Obviously, the edge {u, u’} ¢ F%™” and there is no other path between u and wu’ in the new 
forest. So then, the sum over all possible oriented spanning trees in the left-hand side of 
relation (A.8) will give the weight of all possible forests such that 
Yd wTar(@ > y) A WF) eyez 
y T yY zEÁpry 

where 


Apay = {z E V; A(z Sy) Ee F™} 
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is the set of all vertices from which there is no path to y in the forest F”. Hence, according 


to the first case in Lemma|A.4| the tree-term of V” is 


Vie = = SEE, )qr(z > y). 


yEeEV T 


We are ready for the main result of this Appendix. 


Proposition A.6. The graphical representation of the quasipotential in (A.4) for 
f(y) =P) - (P), Py) =S ky, 2)ay,),  a(a,y) = —a(y, 2) 
consists of two distinct classes of terms, trees and loops: 


V(x) = Vireel £) + Vioop( £) 


where 
Vrela) = > DE T)ar > v) (A-9) 
yEV T 
Vel) E E E O- p DFP. (A10) 
LEL KEKe æ KEO(K) y 


Proof. From Eq. (A.5) we get that 
V(x) = V? (x) — VP (£). 


Now use the representation of V? (x) from Lemma |A.5| The representation of (P) is given 


in Lemma |IV.2 


Lemma A.7. Consider a connected graph G such that the edge (x,x') € E(G) then 

Sow (F>?) => w( FP =y) — 0. 
Proof. Take arbitrary z and F € F”?#. Then F consists of two disconnected trees T} and 
T, where T} is a tree rooted in z, x € Tı and T is a tree rooted in some vertex r. There 
are two possibilities, x’ € T} or x’ € To. If x’ € T, then also F € Fe, If x’ € T, then 


F € F”>". In the same way, if we fix y, every forests F € F”! corresponds to a forests 


F € F*~* for some z. This is a one-to-one correspondence. 
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2. Difference-quasipotential on an edge 


Consider an edge e := {x,y}. We write € := (x,y) when it has a direction. Put V(é) := 
V(x) — V(y). From Proposition (A.6), 


V (2) = Vireel€) + Vioop (2)- 


Recall (A.9) 


Vreel®) = 5 X w(Ta) ar(a > u) —- — SEE Jarly > u) 
uEeV T ueEN T 
= LE where > u) + arlu > y) 
uEeV T 
= y LX oTe sw. 
uEeV T 


According to Proposition[A.6]and Lemma|A.7]the difference of loop terms is given as follows: 


Voow(®) = = E ae M- Lov E we (A.11) 


en KEK æ REO(K LEL KEK ey REO(K) 


where Ky, denotes the set of all tree-loop-trees such that x is located in the tree-loop part. 
If x and y are located in a same tree-loop-tree such that both are in the tree-loop part, then 
the result of equation is zero. We continue with the case that x and y are located 
in a same tree-loop-tree but such that one of them is located on the tree-loop part and the 


other is located in the tree part: 


Vioop (@) = FEl » Di = > ye P0). (A.12) 


LEL KEK ely REO(K KEK, syle REO(K 


Here, Ke .|, denotes a tree-loop-tree with non-oriented loop ¢ such that x is located in the 
tree-loop part and y is located in the tree part. In the case that the edge e is not located on 
a tree-loop Vioop(€) = 0. Put He, as the set of all spanning tree-loops (including the loop £) 
such that the edge e is located on a tree and HO is the set of all tree-loop-trees which are 
created by removing the edge e from a spanning tree-loop of Hee. Corresponding to what 
state of the edge e is closer to the loop, the components of Hee can be split in two groups. 


The set of spanning tree-loops where x is closer to the loop is denoted by He, e, while He 
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is the set of all tree-loop-trees such that x is located on the tree-loop part. Rewrite (A.12) 


as 


Vioop(@) -7(d D a -YS Mw Fall) (A.13) 


LEL Hex?) HeO(H LEL Hen, HeO(H) 


Notice that for every spanning tree-loop including the edge e on a tree either HO or 


A happens. 


As an example, consider the graph in Fig. [14] which has three loops. 


FIG. 14: A graph with three loops. 


The possible tree-loops are shown in Fig. 


b x x 
v u Z y Vv u Z y v u Z y 


(1) (2) (3) 
x x x 
v Uu z y v U z y v u zZz y 
(4) (5) (6) 


FIG. 15: All spanning tree-loops of Fig. h4. 


To find the differences of panne over the edge € the tree-loops including the edge 
e in a tree are engaging, see (2) and (3) in Fig. [15} 
We first look at the unoriented tree-loop Hz. and we remove the edge e. In that way a 


tree-loop-tree is created. Secondly, we consider different possible orientations for the created 
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tree-loop-trees. We rewrite (A.12), 
1 oa 
Wm EE E E nOD) 
£ HE€Hee HeO(H) 


where og(€) = +1. If the edge @ is oriented towards the loop of H, then og(€) = —1; 


otherwise it is positive. So then, 


Vioop(€ S Sooo )w(H), e € trees (A.14) 


W ee. HeEO(H)) 


where He is a set of all spanning tree-loops including the edge e in a tree. If H € He, then 
H denotes a tree-loop-tree which is created by removing the edge e from the tree-loop H. 
O(H?) is the set of all oriented tree-loop-trees made by giving all possible orientations to 
H©) (remember Definition [II.2). Finally, the difference of the quasipotential over a directed 


edge € is 


v(e) = Viree (e) T Vioop (€ (€) 


TED )ar(z > y) 


uEeV T 


+ = > ` on(€)q(@)w(H), e € trees. 


HeHe HeoO(H)) 
Appendix B: Proof of Lemma 


Lemma B.1. For any directed edge € = (£e, Ye), the average of P(x.) = D1, kale) is 


(P= y k(@)a(@)p(ae) = wW S w(A)a(0) (B.1) 
Te Ye HEO(H) 
where w(H) = J [zen K(2, 2’). 
Proof. Take a spanning rooted tree Tz, and the edge € = (ze, y-). We consider two cases: 


First case, if the edge & = (ye, £e) € Tae, then 


k(@)w(Tr.)a(€) = k(@)w(Ty.)a(@) = —k(e')w( Ty. )a(€’). 
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From here we can conclude that 

>, do k@)w(Te.)a@) = 0. 

T €e ETze 
Second case, if the edge @ = (ye, £e) ¢ Tre, then k(€)w(Tz,) is the same as the weight of a 
graphical object made by adding the edge € to the rooted tree Tz.. That graphical object is 


a tree-loop. As a consequence, 


k(@)w(Tx.)q(@) = wale) 


where the loop of H is in the same direction as €. Summing over the edges in loop £ (in 


clockwise or counter clockwise direction) in the underlying graph G gives 


> k@)w(Te,)a(@) = wA). 


ecl 


Finally, we obtain 


Le Ye wW Te Ye 
1 EOIS 
=p È DHAT.) 
= EE E KOOT.) +E E OOUT.) 
E€ Ta? E TZE 
= EANA 
-5 E uA 
HEO(H) 
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